Effects of short-range order on the electronic structure of disordered metallic systems 
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For many years the Korringa-Kohn-Rostoker coherent-potential approximation (KKR-CPA) has 
been widely used to describe the electronic structure of disordered systems based upon a first- 
principles description of the crystal potential. However, as a single-site theory the KKR-CPA is 
unable to account for important environmental effects such as short-range order (SRO) in alloys 
and spin fluctuations in magnets, amongst others. Using the recently devised KKR-NLCPA (where 
NL stands for nonlocal), we show how to remedy this by presenting explicit calculations for the 
effects of SRO on the electronic structure of the bcc CusoZn^o solid solution. 



PACS numbers: 71.15.Ap, 71.23.-k, 71.20.Be, 71. 15. Mb 



I. INTRODUCTION 

Currently, the first-principles theory of electrons in 
disordered metals is based upon Density Functional 
Theory (DFT) and either the Korringa-Kohn-Rostoker 
Coherent-Potential Approximation (KKR-CPA)^* 2 * 3 . or 
its stripped down version, the LMTO-CPA,! for aver- 
aging over the random configurations. This approach 
has been successfully applied to cases where the disor- 
der is internal as well as external. Examples of the lat- 
ter are metallic solid solutions such as Cu c Zri(i- c ) and 
C , M c Pd( 1 „ c ) above their ordering temperatures To. Ex- 
amples of the former are Fe or Ni above their Curie tem- 
peratures T c , where randomness in the crystal potential 
seen by an electron is the consequence of Disordered Lo- 
cal Moments (DLM)^ and secondly valence fluctuating 
systems such as Ce& Despite significant achievements r& 
this methodology suffers from the shortcoming of not 
describing correlations in the fluctuations of the crys- 
tal potential. However a generalization of the KKR- 
CPA theory has recently been proposed, the KKR- 
NLCPA2ii£ (where NL stands for nonlocal), which sys- 
tematically takes into account such correlations, enabling 
environmental effects such as short-range order (SRO) 
to be taken into account. Although the full 3D KKR- 
NLCPA formalism was given in Refs. l9llTol in this pa- 
per we present the first realistic 3D implementation of 
the theorjiii by illustrating the effects of SRO on the 
CusoZn^o system. 

The physics of the above SRO plays a particularly im- 
portant role near phase transitions where it is frequently 
a precursor for long range order and can be said to be 
driving the ordering process. For example, in the order- 
ing of the Cu5oZn5o solid solution into an intermetallic 
compound of B2 symmetry the system lowers its free en- 
ergy by having unlike neighbours more frequently than 
like neighbours even in the disordered state, thereby low- 
ering the temperature To where the system must finally 
order. Such SRO is also central to the understanding of 
electronic transport in general and in iC-state alloys in 
particulariiSii 3 . Moreover, the formation of the moment 



in the DLM state of and the creation of 7-like Ce 
atoms near the 7 — a transition* will be materially af- 
fected by the SRO. The KKR-NLCPA method to be il- 
lustrated here will enable these important problems to be 
tackled in a parameter-independent and material-specific 
way. However, before getting on with the task we com- 
ment briefly on efforts addressing the same problem as 
we do by adopting alternative strategies. 

There have been several attempts to develop cluster 
generalisations of the single-site CPA. As with the CPA, 
the main construct is an effective medium so that the 
motion of an electron through it approximates the mo- 
tion, on the average, of the electron in the disordered sys- 
tem. An early example, the Molecular CPA (MCPA), 1 ! 
introduces a supercell so that the medium has the un- 
satisfactory attribute of broken translational symmetry. 
The Embedded Cluster Method (ECM) 1 ^ refers to the 
non-self consistent embedding of a cluster with all the rel- 
evant disorder configurations into the CPA medium. The 
Travelling Cluster Approximation (TCA)i£ based on di- 
agrammatic methods and the Cluster-CPA (C-CPA)^^ 
based on the Augmented Space Formalism (ASF)22i2i 
are, like the NLCPA, satisfactory on account of their 
translationally invariant, self-consistently determined ef- 
fective media and herglotz analytic properties. Both 
the TCA and C-CPA become rapidly computationally 
intractable, however, and a KKR version of the latter 
has been applied only to model systems ^2*2^ A reason- 
ably good alternative starting point for the electronic 
structure of some disordered alloys is the tight-binding 
(TB)-LMTO methodi combined with the CPA, which 
can include an approximate treatment of the charge self- 
consistency needed for a DFT. Mookerjee and Prasad 23 
have developed a generalised ASF with correlated vari- 
ables to describe SRO, which has been combined with 
the TB-LMTO and real space recursion technique. 24 
This approach has been used successfully to describe ef- 
fects of SRO on the densities of states of several alloy 
systems Nonetheless it is desirable to develop a com- 
putationally tractable generalisation of the CPA within 
the KKR method, the KKR-NLCPA, with fewer approx- 



imations and superiority with regards to accuracy and 
reliability over LMTO methods. It will also be amenable 
for incorporation into a full DFT description of disor- 
dered materials with SROiSl 

This paper is organised as follows. In Sec. II we briefly 
summarise the idea of the KKR-NLCPA (for the full 
derivation see Refs . l9l fToh . and in particular we clearly ex- 
plain how to carry out the fundamental 'coarse-graining' 
procedure for general lattices. Our aim is to show how 
current KKR-based computational codes can be straight- 
forwardly adapted to include the KKR-NLCPA with its 
capability of dealing with disordered systems with SRO. 
In Sec. Ill we present results including SRO calculations 
for the CusoZn^o system, and we conclude in Sec. IV. 



II. FORMALISM 

The first step is to define the scattering path matrix 
t_ ij describing the motion of an electron in an effective 
medium, which ideally should be determined so that it 
would describe the average properties of an electron ex- 
actly. It is a quantity which describes the full effects of 
the coherent potential and is given by 



Real Space 



Reciprocal Space 



t G(R lfc ) + 5G(R lfe ) r 



(1) 



Here a circumflex symbol denotes an effective medium 
quantity and an underscore denotes a matrix in angu- 
lar momentum space. In addition to effective local t- 
matrices t and the usual free-space KKR structure con- 
stants G(Ry) which account for the lattice structure, we 

also have effective structure constant corrections 5G(Ry ) 
which take into acccount all nonlocal scattering correla- 
tions due to the disorder configurations (labelled in 
Ref. Eft. Since the effective medium is translationally- 
invariant, the matrix elements r * 3 are also given by the 
Brillouin zone integral 
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dk[t -G(k)-5G(k) 



ik(Ri-Rj 



(2) 

Since it is not feasible to solve the problem exactly, the 
key idea, based upon concepts from the Dynamical Clus- 
ter Approximation^SSiSS is to perform a consistent coarse- 
graining in both real and reciprocal space in order to 
appropriately deal with <5G(Rj 3 ) and <5G(k) respectively. 
The construction for carrying out this coarse-graining, 
which must retain the translational invariance and point- 
group symmetry of the underlying lattice, has been given 
by Jarrell and KrishnamurthyS? for a 2D square lattice 
in connection with a simple tight-binding model Hamilto- 
nian. We have generalised this construction for realistic 
3D body-centered cubic (free), face-centered cubic (fee), 
and simple cubic (sc) lattices2*iflili which we implement 
for the first time here. 

First we summarise the construction for a general lat- 
tice. Technically, the task is to find an appropriate set of 
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FIG. 1: Example (bec lattice): (a) Cross-section of the 
conventional real-space tiling with Wigner-Seitz cells. The 
shaded sites lie out of the page, (b) Cross-section of the con- 
ventional reciprocal-space tiling with Brillouin zones. The 
shaded sites lie out of the page, (c) Cross-section of the real- 
space tiling (dashed lines) with N c — 2. (d) Cross-section of 
the reciprocal lattice points and reciprocal-space tiles (dashed 
lines) of the coarse-grained system, (e) Cross-section of a real- 
space tile (dashed line) for a N c — 2 cluster containing the 
points Ri = (0,0,0) and R 2 = (a/2, a/2, a/2). The shaded 
sites lie out of the page, (f) Cross-section of the correspond- 
ing reciprocal-space tiles (dashed lines) for the N c = 2 clus- 
ter, with Ki = (0, 0, 0) and K2 = (2n/a, 0, 0) at their centers. 
The shaded points lie out of the page and the solid line de- 
notes a cross-section of the first BZ in the (k x , k y ) plane. The 
BZ can be visualised as a cube with a pyramid attached to 
each of the six faces, and the dotted line shows a projection 
of such a pyramid into the k z plane. 



N c real-space cluster sites {I, J} and corresponding set 
of 'cluster momenta' {K„} satisfying the relation 



J_ _iK n (Rj-Rj 



(3) 
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This may be accomplished as follows: 

• Choose a real-space cluster of N c sites which can be 
surrounded by a tile which a) preserves the point- 
group symmetry of the underlying lattice and b) 
can be periodically repeated to fill out all space. 
For N c = 1, the tiles are the conventional Wigner- 
Seitz cells surrounding each lattice point, as shown 
for the bcc lattice in Fig. [Ja). For N c > 1 there 
may only be solutions to the problem for particular 
values of N c for any given lattice. For the bcc lattice 
the next allowed cluster sizes are N c — 2 and N c = 
16, where the tiles are simple cubes of volume a 3 
and (2a) 3 respectively surrounding each cluster (see 
Fig. flic) for N c = 2). 

• Label the sites of the original lattice by the set 
of vectors {R°™ 9 }, and the centers of the coarse- 
graining tiles by the set of vectors {R^ 9 }. 

• Label the reciprocal lattice corresponding to 
{R°" 9 } by the set of vectors {K°" 9 }. Each K°" 9 
is centered in a Brillouin zone Qbz which periodi- 
cally repeats to fill out all of reciprocal space. For 
the bcc lattice, this Brillouin zone will be a fee 
Wigner-Seitz cell of volume 2(2n/a) 3 , as shown in 

Fig.m». 

• Label the reciprocal lattice corresponding to {R^ 9 } 
by the set of vectors {K^ 9 }. Each K^ 9 is centered 
at a reciprocal space tile (corresponding to the re- 
ciprocal space of the real space tile) which again pe- 
riodically repeats to fill out all of reciprocal space. 
For the bcc lattice, {K- 3 } are simple cubic and will 
be centered at simple cubic tiles of volume (27r/a) 3 
or (tt/o) 3 for the N c = 2 (see Fig^d)) or N c = 16 
cases respectively. 

• Observe that {K°" 9 } C {Kf }. Select N c vectors 
from the set {K^ 9 } which lie within Qbz and do not 
differ by an element of {K°™ 9 }. We define these to 
be the set of cluster momenta {K„}. The recipro- 
cal space contained within the N c reciprocal space 
tiles centered at {K„} is completely equivalent to 
that contained within £Ibz by translation through 
reciprocal lattice vectors {K°" 9 }. See Fig. [Iff) f° r 
the bcc example with N c = 2. 

Refer to the table in Ref.^]for the R/ and K„ values for 
sc, bcc and fee lattices obtained using the above method. 
Note that for the fee lattice the next allowed cluster sizes 
are N c = 4 and N c = 32 respectively, where the real 
space tiles are again simple cubes of volume a 3 and (2a) 3 
respectively surrounding each cluster. See Fig. [3 for the 
fee N c = 4 diagram. 

Having carried out this coarse-graining procedure, we 
can now make an appropriate approximation to deter- 
mine the effective medium. In reciprocal space, we ap- 
proximate 6G(k) within each of the N c tiles by the N c 
values {<5G(K„)}, each defined to be the average of 5G_(k) 
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FIG. 2: (a) Cross-section of a real-space tile (dashed line) 
for a N c = 4 cluster on the fee lattice containing the points 
Ri = (0,0,0), R 2 = (a/2,0, a/2), R 3 = (a/2, a/2,0) and 
R4 = (0, a/2, a/2). The shaded sites lie out of the page, 
(b) Cross-section of the corresponding reciprocal-space tiles 
(dashed lines) for the N c = 4 cluster, with Ki = (0,0,0), 
K 2 = (27r/a,0,0), and K 3 = (0,27r/o,0) shown as the T 
point and the two X points. The fourth tile is centered at 
the X point K4 = (0, 0, 2n/a) and is situated out of the page 
vertically above F. Again the shaded points lie out of the 
page and the solid line denotes a cross-section of the first BZ 
in the (k x ,k y ) plane. 



over the tile centered at K n . The scattering path matrix 
may then be represented by the set of coarse-grained val- 
ues 



l(Kn) 



bz Jn 



dkU' 1 -G(k) -6G{Kr. 



which are straightforward to calculate owing to 6G(J£„) 
being constant within each tile J7k„ • Note that the N c 
integrals here have the same computational cost as one 
standard BZ integral. This is unlike a supercell method 
such as the MCPA where the size of the KKR-matrix, 
which must be inverted at every k-point, increases as N c 
increases. In fact it is straightforward to show that this 
integration step is N c times faster for the KKR-NLCPA 
than for a supercell method for comparable cluster sizes. 
Using Eq. |3J the scattering path matrix at the cluster 
sites becomes 




dkU- 1 -G(k) -SG{K n ) 



xe 



u . 



;K„(R f -Rj) /g\ 

From Nyquist's sampling theorem the effect of 
coarse-graining the effective structure constant correc- 
tions is to reduce their range in real space. In fact from 
Eq. |3| we have 



5G(Ru) = ^^<$G(Kn)e lK " (R '- R '' 
SG(K n ) 



5]£G(R/j)e~ lK " (Rj - Rj) . (6) 
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FIG. 3: (Color online) (a) Density of states (DOS) for pure 
Cu and pure Zn. (b) DOS for ordered Cu^Zn^. 



Note that SG_(Rij) remains a translationally-invariant 
quantity which depends only on the distance between 
sites / and J, now within the range of the cluster size, 
but independent of which site in the lattice is chosen to 
be site I~ It is now straightforward to generalise the 
CPA argument and determine the medium by mapping 
to an impurity cluster problem. We choose a real-space 
cluster consistent with the requirements outlined above, 
and use the Embedded Cluster MethocUSJ£ to replace 
it with an 'impurity' cluster of real t-matrices and free- 
space structure constants in the (still undetermined) ef- 
fective medium. We then consider all paths starting and 
ending on the impurity cluster sites and demand that 
the average over the 2 Nc possible impurity cluster con- 
figurations 7 be equal to the path matrix for the effective 
medium itself i.e. 



7 



(7) 



where -P(7) is the probability of configuration 7 occur- 
ing. Therefore the effective medium t-matrices and ef- 



fective structure constants are determined from a self- 
consistent solution of Eqs. El and EJ An example al- 
gorithm is given in Refs. l9Hlfl SRO may be included 
by appropriately weighting the configurations in Eq. 
(the number of which can be reduced using symmetry 
and sampling), provided that translational-invariance is 
preserved. Observable quantities such as the density of 
states (DOS) can be calculated from the corresponding 
configurationally-averaged Green's function^ The for- 
mula for the DOS within the KKR-NLCPA is given in 
Refs. I9IIT0I Importantly owing to the translational in- 
variance of the KKR-NLCPA medium it is independent 
of lattice site chosen. This is crucial for calculating the 
partially-averaged charge densities to be used in combi- 
nation with DFTi 9 ! 27 ! 32 ! 33 ! 34 Finally, note that the KKR- 
NLCPA formalism reduces to the KKR-CPA for N c = 1, 
and nonlocal scattering correlations (and SRO if desired) 
are systematically included into the effective medium as 
N c is increased, becoming exact as N c — > 00. 



III. RESULTS 

To illustrate the KKR-NLCPA, we present calcula- 
tions for the bec Cu^Zn^ solid solution with lattice con- 
stant 2.86A. In all calculations that follow, Cu and Zn 
potentials come from self-consistent field (SCF) KKR- 
CPA calculations) 3 ^ 3 ^ the Brillouin zone integrals use the 
adaptive quadrature method^ and the energy contour 
has 1 mRy imaginary part. 

Fig. Eta) shows DOS plots for pure Cu and pure 
Zn, and Fig. E{b) shows a calculation for the ordered 
CusoZn^o compound. Since the energies of the Cu and 
Zn d-bands are very different, the system is said to be in 
the 'split band' regime. Physically, this means an elec- 
tron travels more easily between Cu or between Zn sites 
than between unlike sites and so the decrease in over- 
lap between like sites in the ordered case results in a 
narrowing of the Cu and Zn bands by a factor of two 
compared with the pure calculations^ Fig. Eta) shows 
KKR-CPA results for disordered Cu^Zn^Q. It is clear 
that the bands are widened and smoothened compared 
with the DOS for the ordered calculation. The compo- 
nent contributions from Cu and Zn impurity sites em- 
bedded in the KKR-CPA medium are also shown. Next, 
a KKR-NLCPA calculation for disordered Cu^Zn^o for 
a two-site cluster (N c — 2) is shown in Fig. EJb). First 
note that there is little observable difference in the total 
DOS compared with the KKR-CPA calculation. This is 
due to the small size of the cluster, and the difference 
due to the nonlocal scattering correlations shows up in 
detail only on a scale of ±1 state/atom/Ry, as shown in 
Fig. E(c). As expected, it is clear that the extra structure 
is in the energy regions of the impurity d-bands. However 
the most striking aspect of the KKR-NLCPA calculation 
is that the component contributions to the total DOS 
from the four possible cluster configurations are appar- 
ent. The component plots here are the DOS measured 



5 




at the first cluster site when a particular cluster configu- 
ration is embedded in the KKR-NLCPA medium, which 
is the Cu site for the Cu-Cu and Cu-Zn configurations, 
and the Zn site for the Zn-Cu and Zn-Zn configura- 
tions. Crucially, owing to the translational invariance of 
the KKR-NLCPA medium, measurement at the second 
site gives the same results with a simple reversal of the 
labels of the Cu-Zn and Zn-Cu components. These com- 
ponent plots are particularly useful for interpreting the 
effects of SRO on the electronic structure, as described 
in the next section. Finally, Fig. HJb) also shows total 
DOS results for the larger N c — 16 cluster, where the 
extra structure is more visible in the energy region be- 
tween 0.10 Ry and 0.15 Ry where some states are shifted 
to higher energies. 
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FIG. 4: (Color online) (a) Total average density of states 
(DOS) for disordered bcc Cu 5 aZn 5 o using the KKR-CPA. 
Also shown are the contributions from the Cu and Zn compo- 
nents (site restricted average DOS). Ef is the Fermi energy, 
(b) Total average DOS for bcc Cu$oZn$a using the KKR- 
NLCPA with N c — 2, along with the contributions from the 
4 possible cluster configurations (cluster restricted average 
DOS) measured at the first site i.e. Cu for Cu-Cu, Cu-Zn, 
and Zn for Zn-Cu, Zn-Zn. (Owing to the translational in- 
variance, contributions measured at the second site would give 
the same results with a simple reversal of the labels). Also 
shown are total DOS results for N c = 16. (c) Plot of the 
difference in the total DOS between the N c = 1 and N c — 2 
calculations i.e. (total N c = 1) — (total 7V C = 2). 



Short-Range Order 

Unlike single-site theories, it is possible to include the 
effects of short-range order (SRO) on the electronic struc- 
ture of disordered systems using the KKR-NLCPA. This 
may be done by using an appropriate non-random proba- 
bility distribution when averaging over the impurity clus- 
ter configurations in Eq. [7| Of course as the size of the 
cluster increases, the range of possible SRO that may 
be included also increases. Here we show results for bcc 
CusoZn^o using a pair cluster (N c = 2), where it is possi- 
ble to include SRO between nearest neighbour sites only. 
This may be straightforwardly done by introducing the 
nearest neighbour Warren-Cowley SRO parameter 
and using probabilities defined as 



P{CuCu) 
P{ZnZn) 
P(CuZn) 
P(ZnCu) 



P{Cuf + a/A 
P(Zn) 2 + a/A 
P{Cu)P{Zn)-a/A 
P(Zn)P{Cu)-a/4: 



For Cu 50 Zn 50 , P{Cu) = P(Zn) = 0.5 and so the SRO 
parameter can take values in the range — 1 < a < 1, 
where —1, and +1 correspond to ideal ordering, com- 
plete randomness, and ideal clustering. 

Fig- 111 shows the effects of short-range clustering upon 
the configurationally-averaged DOS. It is clear that as a 
increases above zero, the probability of like pair compo- 
nents increases whilst that of unlike pairs decreases, re- 
sulting in corresponding changes to the component and 
total DOS. Indeed at a = +1 the probability of unlike 
pairs is zero and the total DOS is now completely domi- 
nated by the features of the Cu-Cu and Zn-Zn compo- 
nents, as shown in Fig.^f). As expected these features 
are reminiscent of the pure bands shown in Fig.^a), for 
example the magnifying of the trough and peak cither 
side of 0.15 Ry and 0.5 Ry can all be associated with 
the DOS for pure Zn and pure Cu respectively. More- 
over, a new peak appears just before 0.4 Ry associated 
with pure Cu. Fig. [5] shows the effects of short-range 
ordering upon the configurationally-averaged DOS. Here 
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FIG. 5: (Color online) (a) Total and cluster component density of states for bcc Cu^oZn^o using the KKR-NLCPA with N c 
and q = 0. (b)— »(f) Same as (a) but with increasing values of a, corresponding to short-range clustering. 



as a decreases below zero, the components of the total 
DOS due to like pairs decreases whilst that due to unlike 
pairs increases. Finally at a = — 1, as shown in Fig. EJf ) , 
there are only contributions remaining from unlike pairs. 
Evidently, the DOS in general has a closer resemblance 
to that of ordered Cu$oZn$a shown in Fig. |3Jb). For 



example the peaks either side of 0.15 Ry and 0.5 Ry are 
now of roughly equal magnitude and there is also a slight 
overall narrowing of the bands compared with those for 
positive values of a due to the decrease in probability of 
like neighbours. 
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FIG. 6: (Color online) (a) Total and cluster component density of states for bcc Cu^oZn^o using the KKR-NLCPA with N c 
and q = 0. (b)— >(f) Same as (a) but with decreasing values of a, corresponding to short-range ordering. 



IV. CONCLUSIONS 



In summary, the recently devised KKR-NLCPA gen- 
eralises the widely used KKR-CPA method by includ- 
ing SRO, and satisfies all of the requirements for a 



satisfactory effective medium cluster theory^ It has 
relatively low computational cost in comparison with 
supercell-based methods since the BZ integration, the 
most computationally-demanding aspect of a band struc- 
ture calculation, does not scale with the cluster size. In 
this paper we have implemented the KKR-NLCPA for a 



8 



realistic system by illustrating the dramatic changes that 
can occur in the DOS for bcc Cu 5Q Zn 50 . 

Ultimately the KKR-NLCPA will be the electronic 
structure component of a fully self-consistent theory of 
disordered systems. The next step will be to com- 
bine it with DFT, 27 which will also enable charge 
correlations22i2&2i and local lattice displacements^ to 
be systematically taken into account for alloys. Then, 
the SRO parameter a will need to be coarse-grained via 
Eqs. El and could be determined via a linear response 
calculation 41 before being fed back into the electronic 
structure, resulting in a completely ab-initio theory of 



SRO at a given temperature T. Similar treatments will 
also be available for effects such as magnetic SRO in 
metallic magnets at finite temperature. 
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